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. Observations of H2 line emission have revealed higher-than-expected gas temperatures in a 

{^JQ' number of photodissociation fronts. We discuss the heating and cooling processes in photodis- 

i~{ , sociation regions. Observations of NGC 2023 are compared to a theoretical model in which 

• there is substantial gas at temperatures T = 500 — lOOOK heated by photoelectric emission and 

' collisional deexcitation of H2 . In general the model successfully reproduces the observed H2 line 

emission from a wide range of energy levels. The observed [Sill]34.8/im emission appears to 
indicate substantial depletion of Si in NGC 2023. 
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1. Introduction 

00 ' A significant fraction of the ultraviolet radiation emitted by massive stars impinges on 

, the molecular gas associated with star formation. The resulting photodissociation regions 

(PDRs) therefore play an important role in reprocessing the energy flow in star-forming 
galaxies. Modelling these PDRs is therefore an important theoretical challenge, both 
to test our understanding of the physical processes in interstellar gas, and to interpret 
observations of star-forming galaxies. 
O ' It is frequently the case that the illuminating star is hot enough to produce an H II 

^ , region, in which case the photodissociation region is bounded on one side by an ionization 

^ ' front, and on the other by cold molecular gas which has not yet been appreciably affected 

by ultraviolet radiation. The hiy < 13.6eV photons propagating beyond the ionization 
front raise the fractional ionization, photoexcite and photodissociate the H2, and heat 
, the gas via photoemission from dust and collisional deexcitation of vibrationally-excited 

\^ ' H2 . Figure shows the different layers within a PDR. 

2. PDR Thermometry 

The lower rotational levels of H2 tend to be in LTE (except perhaps for the ortho/para 
ratio, which adjusts relatively slowly), so we can use the relative level populations of H2 
as a PDR thermometer. 

We can also use the absolute column densities N{v,J) to tell us how much molecular 
gas is at different temperatures. Because the H2 quadrupole vibration-rotation lines are 
generally optically thin, only dust extinction affects the radiative transfer. 

The level populations N{v, J) can be studied using vibrational transitions (ground- 
based K,II,J,I,R band spectroscopy) or ?;=0-0 pure rotation hues (some from the ground, 
but mainly from space, with ISO or SIRTF, or from the stratosphere, with SOFIA). 



3. Modeling PDRs 

The central physical process in a PDR is the photoexcitation and photodissociation 
of H2 through the Lyman and Werner band lines of H2, and it is therefore important 
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Figure 1. Schematic diagram showing the diflferent zones in a photodissociation region. 



to model this process accurately. In our models (Draine & Bertoldi 1996) we explicitly 
include radiative transfer in the 28765 permitted Lyman and Werner band lines between 
H2 levels with J < 29, using wavelengths and oscillator strengths from Abgrall & Roueff 
(1989) and Abgrall ct al. (1993a, b). The lines are treated independently assuming Voigt 
line profiles and attenuation by dust, but we include a a statistical correction for the 
effects of line overlap of the H2 lines (Draine & Bertoldi 1996). 

Lyman and Werner band photons are absorbed by H2 in the ground electronic state 
X^E+, resulting in excitation to vibration-rotation levels of either the B^S+ or C^ll^ 
electronic states. The electronically-excited state will decay by spontaneous emission 
of an ultraviolet photon. About 85% of the downward transitions will be to a bound 
{v, J) level of the ground electronic state, but about 15% of the transitions will be to the 
vibrational continuum of the ground electronic state, producing two H atoms, typically 
moving apart with a kinetic energy of a fraction of '^l eV (Stephens & Dalgarno 1973). 

Bertoldi & Draine (1996) discussed the propagation of coupled ionization-dissociation 
fronts, and showed that except when the PDR is driven by radiation from a very hot 
star, it is usually the case that the PDR moves into the molecular cloud slowly enough 
that thermal and chemical conditions are close to being in steady-state balance: most 
importantly, H2 destruction is nearly balanced by H2 formation, and heating is nearly 
balanced by cooling. It is th(ui convenient to approximate the structure of the PDR by 
requiring precise thermochemical steady state at each point. The H2 abundance is there- 
fore determined by a balance between destruction by photodissociation and formation 
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Figure 2. N{v,J)/gj for vibrationally-excited levels of H2 toward the southern bar in NGC 
2023. Open symbols: plane-parallel model with nn = 5 x 10*cm~^ and x ~ 5000 observed from 
angle — 78° rol. to surface normal. Filled symbols: level populations obtained from observed 
line intensities after correction for foreground extinction with Ak = 0.5mag. Data from Burton 
et al. (1998) and McCartney et al. (1999). 



on grains, 

-RgrnHn(H) = Cpd"(H2) (3.1) 

where Cprf is the local rate for photodissociation of H2, and Rgr is the rate coefficient 
for H2 formation on grains, with the grain abundance assumed proportional to rifj = 
2n(H2) + n(H) + n(H+). We solve for the steady-state populations of the 299 bound 
(v, J) levels with J < 29, including UV pumping, spontaneous decay, and collisional 
transitions. The distribution over {v, J) of newly- formed H2 produced on grain surfaces 
is uncertain; our models allow us to consider various possibilities. 

The gas is heated primarily by photoelectrons ejected from dust grains, and (for den- 
sities riH ^ lO'^cm"'^) collisional dccxcitation of vibrationally-cxcitcd H2 resulting from 
UV pumping. In our models we use photoelectric heating rates estimated by Weingart- 
ner & Draine (2000a); we assume dust with Ry = 5.5 with 5% of the solar carbon in 
ultrasmall carbonaceous grains. 

The cooling of the gas is dominated by [CII]158/zm, [SiII]35/um, [OI]63/xm, and [FeII]26/um. 
The [OI]63/im line is often optically thick; we use an escape probability treatment. 
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Figure 3. Same as Fig. g, but for = 5 — 8. 
NGC 2023: 78S.11W of HD37903 
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Figure 4. Same as Fig. 0, but for i; = 8 — 12. 
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4. NGC 2023: A Test Case 

The B1.5V star HD 37903 is situated just outside the surface of the L1630 molecular 
cloud, resulting in the famous reflection nebula NGC 2023. Many emission lines of H2 
have now been measured for NGC 2023, making it an ideal test for PDR models. High 
resolution images of NGC 2023 in the l-OS(l) line have recently been published (Rouan 
et al. 1997; Field et al. 1998; McCartney et al. 1999), showing pronounced filamentary 
structure, notably strong emission along the "southern bar" 78" south of HD 37903. 

Black & van Dishoeck (1987) produced the first models to explain the H2 fluorescence 
from NGC 2023. Their favored model assumed cold H2 with a density nn ~ lO^'cm"'^ and 
an ultraviolet radiation fleld enhanced by a factor x ~ 700 relative to the Habing (1968) 
intensity. Draine & Bertoldi (1996) argued that some warm gas with T « 500 — lOOOK 
with density nn ~ lO^cm^"^ was required to explain the observed line ratios. 

New observations have become available both from the ground (Burton et al. 1998; 
McCartney et al. 1999) and from ISO (Bertoldi et al. 2000a). From ground-based spec- 
trophotometry of the bright "southern emission bar" it is now clear that the foreground 
extinction at K is Ak ~ 0.5mag, implying a far-red extinction ^o.S/jm ~ 2.7mag. For 
each observed H2 emission line we compute the column density N(v, J) in the excited 
state required to reproduce the observed surface brightness, for an assumed foreground 
extinction Ak — 0.5mag. 

A detailed discussion of our modelling of NGC 2023 will be reported elsewhere, but 
here we show one example of a model developed to reproduce the observations of the 
southern emission bar. We assume the incident radiation field to be enhanced by a factor 
X ~ 5000 relative to the Habing (1968) fiux, and a gas density nn = 5 x lO'^cm"'^. The 
PDR is assumed to be plane-parallel, inclined relative to the line of sight with cos6' = 0.2, 
where 9 is the angle between the normal to the PDR (= the direction from the PDR to 
the illuminating star) and the direction to the Earth. 

In Figures ||-^ we show N{v, J)/ gj for the model together with the observed level 
populations, versus the energy E{v, J) of the excited state. Open symbols are model 
results; filled symbols are observations. The model is seen to be in quite good agree- 
ment with the populations of the vibrationally-excited levels all the way up to ti = 12! 
The generally excellent agreement between model and observations indicates that the 
description of the ultraviolet pumping process is basically sound. The fact that for a 
given vibrational level the rotational level populations appear to be in good agreement 
with observations indicates that the model has approximately the correct temperature 
and density. Only for a few weak lines are the reported fluxes very different from the 
model predictions: 12-6 Q(5) A = 0.8225^m, where the reported flux is a factor of 20 
stronger than expected, and 9-4 S(ll) A — 0.7663/.tm, where the observed flux is a factor 
~7 stronger than expected. With so many other lines in good agreement, perhaps these 
lines have been misidentified. Note, however, that these are the only two lines from levels 
with E{v, J)/k > 48 X lO'^K, so perhaps there are processes populating these very high 
energy levels which have not been included in the model (see Bertoldi et al. 2000b). 

Figure |5| shows the v = level populations for our model, together with the level 
populations inferred from the ISO observations toward NGC 2023 (Bertoldi et al. 2000a), 
except that the surface brightnesses measured by the relatively large ISO beam have been 
arbitrarily increased by a factor 1.8 - the reasoning here is that the bright southern bar 
represents a region with unusually high limb-brightening, probably about a factor of 2 
higher than the average over the ISO beam. With this adjustment, we obtain quite good 
agreement for most of the rotational levels; the largest discrepancy is for J — 3, where 
the model surface brightness is about twice the observed value. Note that the observed 
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Figure 5. Level populations for rotationally-excited levels of the ground vibrational state in 
NGC 2023. Beam-averaged column densities observed by ISO have been multiplied by 1.8 to 
allow for likely enhancement of the surface brightness on the southern bar. 



emission from the J = 15 level is in excellent agreement with the model! The level 
populations obviously depart strongly from a single-temperature fit; this is mainly due 
to the range of gas temperatures present (see Figure ||) and the excitation of newly- formed 
H2 leaving grain surfaces. 

The model is shown in Figure Photoelectric emission from dust grains dominates the 
heating over much of the PDR, but for 10^^ ^ iVn ^ lO^^cm"^ the dominant heat source 
is collisional deexcitation of H2 which has been vibrationally excited by UV pumping. 
The combination of photoelectric heating and heating by deexcitation of H2 manages to 
maintain a gas temperature of ~ lO'^K over a substantial zone, with the gas temperature 
still as high as ~ TOOK when the H2 fraction has risen to 10%. 

The model assumes gas phase abundances C/H=1.4 x 10^** (1/3 solar), 0/H= 3.2 x 
10-4 (1/2 solar), Si/H=9.0 x 10"^ (1/40 solar), Fe/H=l. 3 x 10"^ (1/250 solar). Predicted 
fine structure line intensities for the bright southern bar (for limb-brightening 1 / cos 6 = 
5) are given in Table |l|. Most of the lines are not yet observed, but the intensity measured 
by ISO for [SiII]34.8/im is only 1/25 of the model prediction! Part of the discrepancy 
can be attributed to beam dilution: the ISO beam is 660arcsec^, whereas the bright bar 
is only ~50arcsec^. However, there presumably is [Si II] emission away from the bar, 
so it appears likely that the Si is depleted by ~400 relative to solar. It is interesting 
to compare this with the S140 PDR, where the inferred Si abundance was 1/40 solar 
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Figure 6. Lower panel: temperature T, ionization fraction rie/riH, and H2 fraction as function 
of column density A^h measured from the ionization front. Upper panel: Heating or cooling 
contributions due to photoelectric emission from dust (p.e.), fine structure emission (f.s.), optical 
line emission (opt), and coUisional oxcitation/dccxcitation of H2 (H2). Solid lines indicate 
heating; broken lines indicate cooling. H2 heats the gas over the region 5000 ^ T ^ 300K. 

(Timmcrmann ct al. 1996). Walmslcy ct al. (1999) have recently discussed the puzzling 
variations in Si abundance observed in PDRs. 

5. Discussion 

There are a number of areas for improvement in our modelling of PDRs: 

• Wc need accurate cross sections for collisional excitation and dcxcitation of H2. 

• We need to understand the details of the H2 formation process in PDRs - the overall 
rate and the {v, J) distribution of newly-formed H2. 

• The rate of photoelectric heating by dust is critical, particularly since drift of the 
dust grains can alter the dust /gas ratio in the PDR (Weingartner & Draine 1999, 2000b). 

• Finally, we need more accurate spectroscophotometry, with refined angular reso- 
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line predicted observed 

[SiII]34.81/ini 1.5 x lO"'^ 5.7 x 10"^ erg cm"^ s"^ 

[FeII]25.99Mm 1.8 x 10"^ 

[OI]63.18pm 4.5 x 10"^ 

[OI]145.5/im 2.9 x 10"^ 

[CII]157.7Atm 4.6 x 10"^ 

Table 1. NGC 2023 Fine Structure Line Emission 



lution, to allow us to understand the geometry of the filaments and sheets which are 
apparent in high-resolution images of NGC 2023 and other PDRs. 

This research was supported in part by NSF grant AST96- 19429. 
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